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Abstract 

We discuss how the coined quantum walk on the line or on the circle 
can be implemented using optical waves. We propose several optical cavity 
configurations for these implementations. 

PACS: 03.67.Lx, 05.40.Fb 

1 Introduction 

Classical implementations of quantum computing algorithms are interesting for 
both practical and fundamental reasons. They allow actual implementations 
that may serve as model systems for their quantum counterparts, and they 
may shed some light on the profound differences between classical and quantum 
waves. 

The concept behind these classical implementations is simple: waves share 
with quantum systems the possibility of creating superposition states p^. Of 
course, the entanglement of separated systems that is specific to quantum me- 
chanics is not present, but such entanglement does not seem to play a central 
role in at least some of the proposed quantum computing algorithms. This 
subject has been considered by a number of authors PP-01- In particular, im- 
plementations of Grover's quantum search algorithm have been proposed [3], 
and actually implemented 0], using only classical resources. It must be pointed 
out however that the lack of entanglement in the classical case manifests through 
the exponential increase of the ammount of classical resources needed with the 
increase of the problem size [1], and it would be more appropriate to speak 
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of simulations rather than actual implementations, we shall come back to this 
point later. 

The so-called quantum random walk (QW in the following) is a process that 
has received a lot of attention in the quantum computing community [5]-|2(J|; see 
the recent review by Julia Kempe [5T] . The QW was introduced by Aharonov et 
al. |S] ten years ago as the quantum analog of the classical random walk (RW for 
short). The attention devoted to QWs is mainly motivated by the fact that as 
some problems can be solved efficiently in classical computation with algorithms 
based on RWs, there is a reasonable hope that QWs could give rise to algorithms 
faster than classical for some specific problems. In work along these lines, Kempe 
[TU| showed that the hitting time of the discrete QW from one corner of an n- 
bit hypercube to the opposite corner is polynomial in the number of steps, n, 
whilst it is exponential in n in the classical case. Subsequently Shenvi et al. ^3] 
showed that a QW can perform the same tasks as Grover's search algorithm, 
and Childs et al. JH] have introduced an algorithm for crossing a special type 
of graph exponentially faster than can be done clasically. Moreover, as classical 
RWs are used in proofs of complexity to place bounds on the difficulty of certain 
problems, then one might hope for similar proofs for quantum algorithms using 

QWs H21- 

QWs can be performed by a quantum computer, and a number of possible 
implementations have been proposed: Travaglione et al. proposed its im- 
plementation in an ion trap quantum computer, Dur et al. |24| consider the 
QW in an optical lattice, and Sanders et al. [35] have suggested its experimen- 
tal realization in cavity quantum electrodynamics. All these proposals consider 
quantum systems. A different approach was followed by Zhao et al. |2fi| . and 
subsequently by Jeong et al. |2Z], who discussed the implementation of QWs 
using single photons and linear optical elements. But it turns out that the 
quantum nature of light does not play any role in the developement of the walk. 
That is, QWs can be performed within classical optics. This has been noticed 
by Hillery et al. wn0 have developed an interferometric analogy of QWs 
and made a suggestion similar to that of Refs. |26H27| . The problem with these 
last proposals is that the number of optical elements grows quickly with the 
number of steps in the QW. 

In a recent paper we have shown that QWs can be understood as an 
interference phenomenon (and thus as a classical phenomenon in the sense that 
quantization does not play any role in its understanding) and can be clasically 
implemented in optical cavities. The advantage of the optical cavity approach 
is that the resources needed (in the sense of the number of optical elements) 
do not depend on the size of the implemented walk. Moreover, we argued 
that the coined QW on the line had indeed been optically implemented by 
Bouwmeester et al. 29 in the context of an optical implementation of the Galton 
board, although the authors did not explicitly make this link. 

In the present paper we extend the results advanced in |28| by considering 
new cavity configurations, by proposing schemes for the implementation of the 
QW on the circle, and by discussing in detail up to what extent the experiment 
of Bouwmeester et al. (2H| constitutes an actual implementaion of the QW. 
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2 The coined quantum walk 



In order to make this paper as self-contained as possible we briefly review here 
some of the main characteristics of QWs. 

Let us first consider the classical RW on the line. In a common version the 
"walker" - the particle or system performing the RW - takes one step to the 
right or to the left depending on the (random) result of tossing a coin. After n 
steps, the probability of finding the walker at a distance m from the origin is 
given by the binomial distribution which, for large n, is well approximated by 
a Gaussian (dashed line in Fig. 1) whose standard deviation is given by y/n. 

In the coined QW, the role of the coin is played by a qubit, that is, by 
a two-level system. In analogy with its classical counterpart, the quantum 
walker moves to the right or to the left depending on the internal state of the 
qubit. The obvious difference with the classical RW is that the coin-qubit 
can be set in a superposition state, and thus the walker can be described by 
probability amplitudes associated with movement both to the right and to the 
left. After each displacement of the walker, the state of the qubit is set to 
a superposition state by means of a suitable unitary transformation, typically 
(but not necessarily) the Hadamard transformation. This operation plays, in 
a sense, the role of the toss of the coin in the RW. Then a new displacement 
occurs, and so on and so forth. 



Figure 1: Probability distribution for n — 200 for both the classical random walk 
(dashed) and quantum walk (continuous) . The initial conditions chosen for calculating 
the QW were Rq.o — l/v2 and I/o,0 = */v2, see Eas. 1111121 . Notice that the 
quantum P m is null for odd m at odd n. We have represented only nonzero values 
and joined them to guiding the eye. 

The probability distribution in the quantum case is very different from the 
classical one in several respects. It resembles the Airy function, as can be seen in 
Fig.l and we discuss elsewhere the asymptotics of the walk distribution in terms 
of Airy functions [2Hj ■ The standard deviation of the QW probability distribu- 
tion is linear in n; that is, the quantum walker walks quadratically faster than 
the classical walker. With increasing time steps, the probability distribution 
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is null for alternatively even and odd values of m. Except for these nulls, the 
probability distribution is quite flat and uniform in its central region. Finally 
we point out that the form of the probability distribution depends strongly on 
the initial state of the coin. Several analytical approaches have been developed 
that provide explicit solutions of the QW on the line [51151119). 

The QW on the circle is similar to the QW on the line. But the circle 
contains only a finite number of discrete points from m = to m — M — 1 and 
then position m = M + 1 is the same as position and so on. The probability 
distribution is, of course, very different to that of the QW on the line (it does 
not converge to a limit), and in Fig. 2 we have plotted three examples that are 
detailed in the caption (see also, e.g., [9"lll7|'). 



Figure 2: Probability distribution for the indicated number of steps, n, for the quan- 
tum walk on a circle with 61 discrete positions. The initial is i?3o,o = 1/V2 and 
^30,0 = i/V%- 

So-called quantum random walks are of course not random, since the evo- 
lution of the quantum walk is completely deterministic, and randomness enters 
only when a position measurement of the walker is finally performed. Thus we 
refer to this process simply as a "quantum walk." We mention that a different 
but related type of quantum walk, the so-called continuous QW, was introduced 
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by Farhi and Gutmann [7|. The continuous QW does not need the use of any 
coin or qubit. We refer the interested reader to the literature Q |^| for the 
details of this walk, which we will not consider here. 



3 Classical implementation of the quantum walk 
in optical cavities 

3.1 The Hadamard walk 

Let us first consider the implementation of the standard Hadamard walk. In 
our approach, the role of the walker will be played by the frequency of a light 
field, and the role of the coin will be played by some other degree of freedom of 
the light field. This can be its spatial path or its polarization state. That is, 
in the terminology of PP, we consider both position and polarization cebits. A 
cebit is a two-component complex vector representing the light field, and is the 
classical counterpart of the qubit. 
For example, the light field 

i 

E=Y^ E m e l[ ^ 0+ma)t ~ kmZ] +c.c, (1) 

m— — I 

where uj is the carrier frequency, Co is the frequency difference between succes- 
sive frequency components, and z is the propagation direction (the vector E m 
lies in the (x,y) plane), can be represented by the abstract state 

i 

IVO = ^2 [Rm\m,x) + L m \m,y)}, (2) 



m=—l 



where R m = x-u m , L m = y-u m (u m = E m j E m ) and J2 m =-i 
= 1. In this notation the basis vectors \m, c) label the frequency and polariza- 
tion, with c = x,y and we can associate x(y) with the coin head (tail). We can 
also consider position (or direction) cebits, that is, a field of the form 

i 

E=J2 e c E m ^ c e l[{u0+mO)t - k ^ ] + c.c, (3) 

c=l,2 m=-l 

where r c — k c -r (k c denotes the direction of propagation of the c th field), and 
e c denotes its polarization. This can be represented by the state 

i 

|V) = ]T [Rm\ m,ri) +L m |m,r 2 )] , (4) 

m— — l 

where now R m = E„ hl / (\E m> i\ 2 + \E m ^ 2 \ 2 ^ andi m = E ma / (j£m,i| 2 + |^m,2| 2 )- 
Generally we use the index c in |m, c) to label the variable that plays the role 



\Rm\ 2 + \L m \ 2 
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of the internal state of the coin; this is the polarization c = x, y in our first 
example, and the different paths, c = r%, in our second. Of course more than 
two optical paths can be considered, and thus coin-cebits with more that two 
components; this would be useful for the implementation of multidimensional 
QWs 1 1 3 j . but we shall not consider these more complication dynamics here. 

As stated, the walk is executed by the light frequency. Then there must be 
a unitary operator that performs the conditional displacement operation 



V \m, Ci) = |m ± 1, Cj) , i = l,2. 



(5) 



For polarization cebits, V can be physically implemented with an electrooptic 
modulator (EOM) to which a linearly time dependent voltage is applied: the x 
(y) polarization component of the field frequency component (uj + nuj) will see 
its frequency increased (decreased) by an ammount u). For a position cebit, the 
EOMs must be adjusted so they are not polarization selective. Two such EOMs 
are needed, one to increase the frequencies in one of the beams, and the second 
to decrease the frequencies in the other beam. 

After each displacement of the field frequency, a Hadamard transformation 
acting on the cebit has to be implemented. The Hadamard transformation is 
represented by the operator 



H 



1 
71 



1 1 

1 -1 



whose action on the cebits reads 



H\m,d) = -j= [\m,ci) 
H\m,c 2 ) = [|to,ci) 



\m,c 2 )} , 
\m,c 2 )} . 



(0 

(7) 
(8) 



For polarization cebits, the Hadamard transformation is optically implemented 
using a half-wave plate (HWP) with its fast axis forming an angle n/8 with 
respect to the x axis. For position cebits, it is implemented using a beam- 
splitter with two 7r/2 phase shifters 

The QW is implemented by the repeated action on the cebit of the operator 
HV, that is, after n iterations the state reads 



IV0„ 



HV 



IVOo: 



+n 



that can be written as 

W) n = X] \ R m\ m ,Cl) + L m \m,C 2 )] 
m=—n 

with 



L't: 



1 

i 

71 



[■fitn-l.Ti-l + Wl,n-l] i 
[■fitn-l.Ti-l — i 



(9) 
(10) 

(11) 
(12) 
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where R m ,o = £ m ,o = 0, if m ^ and R m ,~i = £ ra ,-i = Vm. Finally, the 
intensity of each frequency component of the light field, which is the optical 
analog of the probability of finding the walker at position m at iteration (time) 
n, is given by 

P m (n) = |i? m ,„| 2 + \L m . n \ 2 , (13) 

which is represented in Fig. 1 1 . 

In order to implement n steps the best option is to introduce the described 
elements in an optical cavity. In the next subsection we propose several designs 
for optical cavities in which the QW on the line is performed. 

3.2 Optical cavities for the QW on the line 

An optical cavity imposes the constraint that the optical frequencies must fit 
within a set of cigcnfrcquencics defined e.g. by the length of the cavity. Thus, 
the time dependent electric field applied to the EOMs and the cavity length 
must be adjusted in such a way that the frequency shift u> matches the cavity 
free spectral range. We discuss this important point below. It is also important 
to notice that the field injected in the cavity from the outside (the initial state) 
can be injected only once to effect the kind of initial conditions described above; 
more complicated walks would be generated from multiple injections. In our 
analysis we assume there are negligible cavity losses; in practice such losses 
would limit the number of steps that could be implemented, unless an amplifier 
were added to the cavity to compensate them. In Figs. 3 and 4 we illustrate in 
schematic form several possible designs for the implementation of the QW on 
the line. 

In Figs. 3(a) and 3(b) a polarization cebit is considered. In the scheme of 
Fig. 3 (a) we represent a ring cavity containing an EOM and a HWP, as discussed 
in the previous section. In the scheme of Fig. 3(b) the concept is the same but 
utilizes a linear cavity. As a result, (i) the HWP has to be replaced by a quarter 
wave-plate, QWP, and (ii) the frequency shift introduced by the EOM must be 
half the size as compared to the ring cavity, as light crosses both elements twice 
in every roundtrip. 

Let us discuss in more detail the value of the frequency jump u>. We have 
commented above that the cavity imposses that u> fits the cavity free spectral 
range (FSR), but the frequency shift introduced by the EOM cannot be better 
resolved than the inverse of the cavity roundtrip time r c , which is precisely 
FSR. Then the steps of the QW would not be well resolved if their size equals 
FSR. Nevertheless there are several ways of overcoming this difficulty. One way 
is to take the frequency jump equal to several free spectral ranges, i.e., Cu = f 
FSR with / an integer larger than one. In this way the uncertainty in the 
frequency displacement does not avoid the resolution of the QW steps. Another 
possibility is that the frequency shift introduced by the EOM be smaller than 

x We notice that the QW can also be defined as the repeated action of the operator VH as 
in 0. In this case the probability distribution is the same but the equations of evolution of 
the amplitudes are different although can be easily related to the ones used here. 
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FSR. In this case it would be necessary several roundtrips to perform a single 
frequency step of the QW and then it would be necessary to control the action 
of the HWP (or the QWP) that performs the Hadamard transformation as it 
should not act until the frequency step of the QW is completed. For example, if 
the frequency step takes five cavity roundtrips, the HWP should act only once 
every five rountrips. This can be accomplished by substituting the HWP by a 
second EOM to which a constant voltage of appropriate magnitude is applied, 
this voltage acting only every five cavity roundtrips. There is still another 
possibility by using a single pulse light field with a pulse duration shorter than 
t c . In this case Q should be larger than the pulse spectral width but it would 
be not necessary to match an integer of FSR. In any case, it is clear that 
the sequence of frequency displacement and Hadamard transformation can be 
performed within an optical cavity with enough resolution. 



Figure 3: Schemes for the optical implementation of the QW on the line using po- 
larization cebits. In (a), a ring cavity, the electrooptic modulator (EOM) shifts up 
(down) the frequency of the x (y) polarization component of the light field in ui, and a 
half- wave plate (HWP) with its axis forming an angle ir/8 with respect to the x-axis, 
performs the Hadamard transformation. In (b), a linear cavity, the EOM modula- 
tor shifts the field frequency uj/2 depending on the field frequency as in (a), and 
a quarter-wave plate (QWP) performs the Hadamard transformation (notice that in 
the Fabry-Perot cavity the light passes twice through each intracavity element every 
roundtrip) . 

In the ring cavity shown in Fig. 4(a), a position cebit is used. Two electroop- 
tic modulators, EOM1 and EOM2, are needed here, each one of them acting on 
each of the two sides of the cavity, which play the role of the two coin states. 
The action of these EOMs must be polarization independent. Alternately, the 
polarization of the light must be aligned with the axes of the EOMs, and the set 
of axes of one EOM must be aligned with the set of axes of the other. One of 
the EOMs increases the field frequency by Q, while the other decreases it by the 
same amount. The Hadamard transformation is performed by the beamsplitter 
BS; we do not represent in the figure the additional phase filters necessary for 
making a Hadamard, 

In Fig. 4(b), again a ring cavity, a hybrid between a polarization and a direc- 
tion cebit is considered. In this bidirectional cavity the horizontal and vertical 
polarization components of the field are forced to travel in counterpropagating 
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Figure 4: Schemes for the optical implementation of the QW on the line using position 
cebits. In (a), two coupled unidirectional ring cavities, EOM1 and EOM2 increase 
and decrease, respectively, the field frequency every roundtrip and the Hadamard 
transformation is performed by the beamsplitter BS. The bidirectional ring cavity in 
(b) is designed for sustaining a x (y) polarized field in the clockwise (counterclockwise) 
direction. The EOM increases (decreases) the x (y) polarized field frequency and the 
set formed by the two quarter-wave plates QWP1 (axis 7r/4), QWP2 (axis — 7r/4), 
and BS performs the Hadamard transformation, see text. 

directions. As in the cases described in Fig. 3, a single EOM is needed that in- 
creases (decreases) the field frequency of the x (y) polarization component in ui. 
The set of devices formed by the two quarter-wave plates QWP1 (axis oriented 
at 7r/4) and QWP2 (axis oriented at — 7r/4) and the BS plays a double func- 
tion: (i) it guarantees that the clockwise (counterclockwise) field remains x (y) 
polarized, and (ii) it performs the Hadamard transformation. The Hadamard 
transformation is implemented by the BS; when done in the presence of the 
two QWP, this guarantees that the transmitted field maintains its polarization, 
while the reflected field changes its polarization from x to y and vice versa. Here 
is how it works: Note that when traversing the three elements the polarization 
of a beam does not change, since the effect of QWP1 is cancelled out by the 
effect of QWP2 , and vice versa. For light reflected by the BS, on the other 
hand, QWP1 (or QWP2) is crossed twice, which is equivalent to the effect of 
a HWP with axis at 7r/4 (— 7t/4); this changes x polarization to y polarization 
(or y to x). 



9 



An interesting aspect of the two schemes shown in Fig. 4 is that the output 
corresponds to only one of the two cebit states and thus the spectrum of the light 
does not correspond to P m (n), Eq.|JT3}, but to P^ (n) = |-R m ,n| 2 or P^ (n) = 
|£m,n| • I n or( i er to obtain the complete QW, one should allow the light to 
exit the cavity through two of the cavity mirrors and then combine the two 
beams into a single beam. In Fig. 5 we represent P m (n) together with P^ (n) 
and P^ (n) for an initial condition different to that of Fig.l, see caption. 



Figure 5: Probability distribution for n = 200 for the quantum walks. The initial 
conditions chosen for calculating the QW were i?o,o = 1 and Lo,0 = 0, see Eas. 11111211 . 
Notice that P m is null for odd m at odd n. We have represented only nonzero values, 
r^i P^i and P m are represented in (a), (b), and (c), respectively. 



3.3 More general walks 

The Hadamard walk is only one of several possible ways of performing a QW. 
What is actually needed is a unitary transformation that mixes the two coin- 
cebit states after each system displacement but it is not necessary that this 
unitary transformation be the Hadamard transformation. That is, there is a 
family of unitary transformations that provide the same final probability distri- 
bution, (see, e.g., |3EJE1)- In particular, Konno ^Jj has studied QWs for 
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which the unitary transformation applied to the qubit after each displacement 
reads 

"-OS)- ,14 » 

where a, b, c and d are complex, in general, and satisfy 

|a| 2 + |6| 2 = | C | 2 + |d| 2 = l, 
ac* + bd* = 0, 

|A| = 1, A = ad -be, 
c= -Ab*,d = Aa*. (15) 

Konno JI] demonstrated that, after a large number of iterations n, for QWs 
governed by such a U the mean displacement of the walker, (#), and the mean 
quadratic displacement, (x 2 ), are given by 



(x) 



a 2 2Re(ab*aP* 

\P\ - M +- 



lal 2 



■(l-|6|)n, (16) 
(x 2 ) = (l-\b\)n 2 , (17) 

where |c) = coZ is the initial state of the coin-qubit. The above results 

apply, obviously, to the Hadamard walk by taking a = b = c = —d = l/y/2. 

For the implementations of the QW here considered, this means that it is 
not essential that the unitary transformation after each frequency displacement 
be implemented by a HWP, in the case of polarization cebits, nor by a 50/50 
beamsplitter for position cebits, as they can be substituted by any unitary 
transformation of the type U on the polarization or position. Now we have the 
elements required for discussing the experiemnt of Rcf. 29 . 



3.4 The experiment of Bouwmeester et al. |29| 

Bouwmeester |29| proposed an optical implementation of the Galton board (the 
quincunx), and studied it both theoretically and experimentally. What was 
actually implemented is a grid of Landau-Zener crossings through which a light 
beam propagates, and the focus was the study of the appearance of recurrences 
in the spectrum of the light field. We refer the interested reader to Ref. 29 for 
details. 

A simplified version of the experimental device used in |29j is that repre- 
sented in Fig. 3(b), but replacing the QWP by a second EOM with its axis 
rotated it /4 with respect to the axis of the first EOM. The second EOM intro- 
duces a dephasing 5/2 between the two polarization components, where 5 can 
be varied at will. Then, after a cavity roundtrip, the matrix representing the 
action of the second EOM action on the x and y polarization components of the 
field reads 

U 5 =( C0S ^ " Zsin /V (18) 
V -iBino coso J 
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It is easy to see that Us can be seen as a particular case of U 114|) for a = d = 
cos S and b = c = — ism 5. From Eas. (|16l) and (|17fl we find 



( x ) = |/?| 2 -|a| 2 + 2Im(a/r)taii(5 

• (1 — sin 5) n, 
(x 2 ) = (1 - sin S) n 2 . 



(19) 
(20) 



Notice that for a = (3 both real quantities - that is, for an initial linearly 
polarized state - (x) = for all 5 and the variance of the spectrum is a = 
vl — sin 5. 

From this it is clear that the experimental device used by Bouwmeester 
can implement QWs. In fact, in Fig. 6 of their work which corresponds to 
S = 7r/5, a power spectrum very similar to those expected for the optical imple- 
mentation of the QW can be clearly appreciated. The authors (2H| considered 
Fig. 6 in their manuscript as a demonstration of the coherence quality of their 
system, and did not make an explicit connection of their work with QWs. 

Nevertheless there is subtle but essential difference between the experiement 
of [3!2| and an actual implementation of the QW. In the frequency displace- 
ment w is a fraction of the FSR but the unitary transformation (18) is applied 
at every roundtrip (see our discussion above in Subsection 3.2). Then, the ex- 
periment in |2H| cannot be regarded as an implementation of the QW although 
is something certainly very similar. That is, the optical Galton board is different 
from the QW. Nevertheless we can consider the experiment of Bouwmeester et 
al. [35] as a proof of principle that the QW can be implemented in an optical 
cavity. 

3.5 Quantum walk on the circle 

The QW on the circle, again with the frequency playing the role of position, 
can be also implemented optically. If 2M — 1 is the number of possible discrete 
values for the frequency, it is necessary to devise a way of shifting the frequency 
u> + (M + 1) Q to u> — Muj and the frequency ui — (M + 1) Q to ui + Mu>. This 
can be done by substituting the EOMs in the schemes of Figs. 3 and 4 by the 
device shown in Fig. 6, which we denote as EOM . It consists of two electrooptic 
modulators and two specially designed mirrors. Mirrors Ml and M2 only reflect 
the frequencies u> + Muj and u> — Muj 1 respectively, and are transparent to the 
rest of the frequencies in the spectrum. 

To see how the circle walk results, we consider as an example the scheme 
in Fig. 4(a). First we look at the left side of the cavity, in which frequencies 
were increased at EOM1 in the QW on the line. With EOM of Fig. 6 replacing 
EOM1 in Fig. 4(a), consider any step at which the field entering EOM contains 
all the allowed frequencies from uj — Muj to u> + Muj. The role of Ml is to 
separate the frequency ui + Muj, which is directed to EOMa, from the rest of 
frequencies, which are directed to EOMb. The modulators are configured such 
that EOMa decreases the entering frequency by 2Mw, while EOMb increases 
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Figure 6: Scheme of the EOM device that substitutes the electrooptic modulators 
in Fig. 3(a) in order to perform the QW on the circle. Ml and M2 are mirrors that 
reflect a single frequency (Mui and — Mui, respectively) and EOMa and EOMb are 
electrooptic modulators. See text for a detailed description of the device operation. 

the frequencies entering it in Q. After traversing the electrooptic modulators, 
the frequencies enter M2. The frequencies {ui — (M — 1) u>, . . . , us + Muj} that 
come from EOMb traverse M2, while the frequency us — Muj, that comes from 
EOMa is reflected by M2. Then, the set of frequencies exiting EOM is the same 
set that entered, but all of them have been shifted appropriately. 

In the right side of the cavity in Fig. 4(a), where frequencies were decreased 
in the QW on the line, EOM2 is substituted by another EOM that works in a 
slightly different way to the one in the left side of the cavity: first, the positions 
of M2 and Ml are interchanged (or, what is the same, light enters from below), 
and second, EOMa and EOMb in EOM are modified in the sense that EOMa 
now increases the incoming frequency in 2Mu> whilst EOMb decreases in us the 
rest of frequencies. In this way, in the right side of the cavity all frequencies 
entering EOM are decreased except us — Muj that is converted into us + Muj as 
it must be. Of course it is essential that the optical lengths of the two paths in 
EOM be symmetric in order not to introduce any spurious dephasing during 
the separation of frequencies. 

We finally comment that the device EOM can be also used in the cavities 
of Figs. 3(a) and 4(b) by suitably adapting its operation. In particular, in the 
scheme of Fig. 4(b) a single EOM is sufficient for performing the QW on the cir- 
cle: as the clockwise and counterclockwise fields have orthogonal polarizations, 
the electrooptic modulators in EOM can increase (decrease) the frequencies 
oppositely for the two polarization components. 

4 Discussion 

We have shown in detail how quantum walks, both on the line and on the circle, 
can be implemented in optical cavities. We have also discussed the close relation 
existing between the optical Galton board of Bouwmeester et al. [21] and the 
QW on the line. In fact we have shown that the experiment of j2HJ can be 
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regarded as a proof of principle that our proposals for the implementation of 
the QW can be carried out in optical cavities. All of this leads to the obvious 
conclusion, as noted earlier by us |2E], that the coined QW does not require 
quantization for its implementation and thus, in this sense, it is not quantum. 
Nevertheless this point requires further discussion. 

First of all we should carefully distinguish between an actual implemen- 
tation and a simulation of the QW. A simulation would not capture all the 
characteristics of the QW whilst an implementation would do. Then, in order 
to make an assesment as strong as "the QW is not quantum but classical" we 
must be completely sure that what we are proposing is an implementation and 
not a simulation as in this last case some quantum features would be lost an 
this would manifest in some measurable characteristics. In order to make this 
distinction clearer, let us consider the optical implementations of the Grover's 
search algorithm 00]: it is rather clear that these cannot be regarded as actual 
implementations but only as simulations as the lack of quantum entanglement 
manifests in the exponential increase of classical resources with the increase of 
the problem size. In other words, with the only use of classical resources, quan- 
tum entanglement can be simulated but not implemented and this manifests in 
the scaling of the needed resources. Then, our proposal is an implementation 
or a simulation? 

In order to answer this question we should first clarify the role played by 
entanglement in the QW, something that is beyond the scope of this article and 
that we plan to investigate elsewhere, but even without this knowledge we can 
still clarify to some extent this point. It is obvious that the classical resources 
needed in the schemes here proposed for the implementation of the QW on the 
line (or on the circle) do not grow with the number of steps in the QW except, 
perhaps, in the detection stage but here the resources needed would grow linearly 
with the number of steps, as would grow in a quantum implementation. Then it 
seems that entanglement does not play any role in the QW and we can consider 
our proposals as actual implementations. But this first conclussion could be 
misleading as the increase in the problem size does not lie in the number of steps 
in the walk but in the dimension of the QW. We have noted above that our 
schemes could be generalized for implementing multidimensional QWs for 
that, coins with more than two states are needed or, in other words, several coins 
are needed. Polarization and direction cebits could be used simultaneously for 
implementing these multidimensional QWs. It would be here where the problem 
of the scaling of the classical resources would appear: a d dimensional QW 
requires d coin-qubits but their classical simulation requires 2 d light beams pQ. 
Then, it is in these higher dimensional walks that quantum entanglement would 
manifest. In conclussion: one-dimensional (even two-dimensional) QWs can be 
classically implemented but higher dimensional QWs can only be simulated. Of 
course this is a very preliminary conclussion as in some cases multidimensional 
QWs could be reduced to equivalent unidimensional QWs 

Let us finally note that one of the interesting phenomena that could be 
studied by using any of the designs here proposed is the effect of decoherence on 
the QW. It has been shown by Kendon and Tregenna ^B] that decoherence can 
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be beneficial for QWs in the sense that it can lead, under suitable conditions, 
to more homogeneous probability distributions. Decohcrence can act on the 
walker and/or on the coin, and it could be simulated in the schemes devised 
here by introducing, e.g., random changes in the polarization or direction of the 
light beams. 



5 Acknowledgements 

This work has been supported in part by the UK Engineering and Physical 
Sciences Research Council and the European Union. ER acknowledges a grant 
from the Ministerio de Education, Cultura y Deportes of the Spanish Goverment 
(Grant PR20002-0244) and partial fmantial support from the Spanish Ministerio 
de Ciencia y Tecnologia and European FEDER funds through project BFM2002- 
04369-C04-01. JES acknowledges financial support from the Natural Sciences 
and Engineering Research Council of Canada. We gratefully acknowledge Dr. 
Javier Garcia (Universitat de Valencia, Spain) for his careful reading of the 
manuscript. 



References 



[< 

[5 
[6 
[T 
[8 
[9 

[10; 
[11 



R.J.C. Spreeuw, Phys. Rev. A 63, 062302 (2001) 

N.J. Cerf, C. Adami, and P.G. Kwiat, Phys. Rev. A 57, R1477 (1998) 

P.G. Kwiat, J.R. Mitchell, P.D.D. Schwindt, and A.G. White, J. Mod. Opt. 
47, 157 (2000) 

N. Bhattacharya, H.B. van Linden van den Heuvell, and R.J.C. Spreeuw, 
Phys. Rev. Lett. 88, 137901 (2002) 

Y. Aharonov, L. Davidovich, and N. Zagury, Phys. rev. A 48, 1687 (1993) 

D. Meyer, J. Stat. Phys. 85, 551 (1996) 

E. Farhi and S. Gutman, Phys. Rev. A 58, 915 (1998) 
A. Nayak, and A. Vishwanath, e-print quant-ph/0010117 



D. Aharonov, A. Ambainis, J. Kempe, and U. Vazirani, Proceedings of 
the 30th Anual ACM Symposium on Theory of Computation (ACM Press, 
New York, 2001) 50; e-print |quant-ph/0012090| 

J. Kempe, e-print quant-ph/0205083 

N. Konno, Quantum Information Processing 1, 345 (2002); e-print 



quant-ph/0206103 



[12] E. Bach, S. Coppersmith, M.P. Goldschen, R. Joynt, and J. Watrous, e- 
print |quanT-ph/0207008 



15 



[13] T.D. Mackay, S.D. Barltctt, L.T. Stephenson, and B.C. Sanders, J. Phys. 
A: Math. Gen. 35, 2745 (2002) 



[14 

[15; 
[ie; 

[17 

[is; 

[19 
[20 

[21 

[22; 

[23 
[24 

[25; 

[26 

[27 

[28; 
[29 



N. Shenvi, J. Kempe, and K.B. Whaley, Phys. Rev. A 67, 052307 (2003) 

A.M. Childs, R. Cleve, E. Deotto, E. Farhi, S. Gutman, and D. A. Spiclman, 
|quant-ph/ 0209131 

V. Kcndon and B. Tregenna, Phys. Rev. A 67, 042315 (2003) 

M. Bednarska, A. Grudka, P. Kurzynski, T. Luczak, and A. Wojcik, e-print 
|quant-ph/0304113| 



M. Hillery, J. Bergou, and E. Feldman, e-print quant-ph/0302161 



H.A. Carteret, M.E.H. Ismail, and B. Richmond, e-print quant-ph/0303105 



B. Tregenna, W. Flanagan, R. Maile, and V. Kendon, e-print 
|quant-ph/0 304204 

J. Kempe, Contemp. Phys. 44, 307 (2003) 

We thank Steven Bartlett, of Macquarie University, for pointing this out 
to us. See, for example, C.H. Papadimitriou, Computational Complexity 
(Addison Wesley, 1994) 

B.C. Travaglione and G.J. Milburn, Phys. Rev. A 65, 032310 (2002) 

W. Diir, R. Raussendorf, V.M. Kcndon, and H.-J. Briegel, Phys. Rev. A 

66, 052319 (2002) 

B.C. Sanders, S.D. Bartlett, B. Tregenna, and P.L. Knight, Phys. Rev. A 

67, 042305 (2003) 

Z. Zhao, J. Lu, H. Li, T. Yang, Z.-B Chen, and J.-W. Pan, e-print 
|quant-ph/0212149| 

H. Jeong, M. Paternostro, and M.S. Kim, e-print quant- ph/0305008 

P.L. Knight, E. Roldan, and J.E. Sipe, Phys. Rev. A 68, 020301 (2003) 

D. Bouwmeester, I. Marzoli, G.P. Karman, W. Schleich, and J. P. Woerd- 
man, Phys. Rev. A 61, 013410 (2000) 



1G 




m 



0.10 




(a) 







E0M1 








(b) 







QWP1 



QWP2 





T 



BS 



0.08 



P R m 0.04 W 



0.00 
0.02 




P m 0.01 



0.00 
0.10 




P m 0.05 



0.00 




-100 







100 



m 




EOM 



